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Natural and artificial light harvesting systems often operate in a regime where the flux of photons 
is relatively low. Besides absorbing as many photons as possible it is therefore paramount to prevent 
excitons from annihilation via photon re-emission until they have undergone an irreversible energy 
conversion process. Taking inspiration from photosynthetic antenna structures, we here consider 
ring-like systems and introduce a class of states we call ratchets: excited states capable of absorbing 
but not emitting light. This allows our antennae to absorb further photons whilst retaining the 
excitations from those that have already been captured. Simulations for a ring of four sites reveal a 
peak power enhancement by up to a factor of 35 under ambient conditions owing to a combination 
of ratcheting and the prevention of emission through dark-state population. In the slow extraction 
limit the achievable power enhancement due to ratcheting alone exceeds 20%. 


Introduction - The absorption of light and prevention 
of its reemission is essential for the efficient operation 
of solar energy harvesting devices [T]. From the many 
causes of device inefficiency, few are as fundamental and 
seemingly insurmountable as energy loss via radiative re¬ 
combination: any absorption process must have a com¬ 
panion emission process. This inherent absorption ineffi¬ 
ciency is a result of the principle of ‘detailed balance’ and 
constitutes a key contribution to the famous Shockley- 
Queisser limit 2|. However, pioneering work by Scully 
showed that it is possible to break detailed balance, given 
an external source of coherence fjj; later work showed 
that this can be achieved by clever internal design alone, 
by using an optically dark state to prevent exciton re¬ 
combination m- Such dark states are populated pas¬ 
sively if the energy separation between dark and bright 
states falls into the vibrational spectrum of the absorb¬ 
ing nanostructure: dissipation then preferentially medi¬ 
ates transfer into states from which optical decay cannot 
occur. This is thought to play a role in photosynthetic 
light harvesters |1 El, e.g. by means of dynamic localisa¬ 
tion reducing the effective optical dipole strength nnj. 

However, time spent in dark states is ‘dead time’ with 
respect to absorbing further photons. In organic light 
harvesting systems the time needed to extract an exci¬ 
ton from a dark state and turn it into useful energy is 
often orders of magnitude slower than typical light ab¬ 
sorption rates mm2!, and this results in the loss of any 
subsequently arriving photons. In this Letter, we show 
that certain exciton states, which we will call ratchets , 
can enhance photocell efficiency by enabling absorption 
of these subsequent photons, while preventing emission. 

Model - We consider a ring of N identical two-level 
optical emitters (see Fig. [I]), with nearest pairs coupled 
by a transition dipole-dipole interaction, governed by the 


Hamiltonian (h = 1): 

N N 

H s a t+ s °*+i + CT i+ i a i ) > (!) 

t=1 i = 1 

where ui is the bare transition energy of the sites, S is 
the hopping strength, the af denote the usual raising 
and lowering operators which create and destroy an exci¬ 
tation on site i , and <jn+ i = o\. For a single exciton, the 
eigenstates are equal superposition of excitons localized 
at the different ring positions, with each being character¬ 
ized by a relative phase hi = 2irj/N with j £ 0, 1...N — 1 
between exciton states in adjacent positions. Full diag- 
onalization can be achieved by the successive applica¬ 
tion of the Jordan-Wigner and Fourier transforms IDE], 
whence the Hamiltonian becomes H s = "Yha ^k \K) (K\. 
K describes the eigenstate with energy A# and each 
K now corresponds to a set of the phase factors ki , 
one for each excitation [see Supplementary Information 
(SI) [15]]. The eigenstates fall into a series of bands, 
where each state within a band contains the same num¬ 
ber of excitons n = \K\, with \K\ the number of elements 
in set K. 

Optical transitions connect eigenstates differing by one 
exciton with rates proportional to F k,k' = \( K '\ J+ \K)\ 
where J ± = y~] ■ 1 erf . An explicit analytical expression 
for Y k ,k- can be derived mm and is given in \m- 0nl y 
the fully symmetric ki = 0 single exciton eigenstate has 
a dipole matrix element with the ground state, since here 
the transition dipoles interfere constructively. The other 
N — 1 single exciton eigenstates do not couple to the 
ground state via light with a wavelength that is much 
longer than the separation between sites. However, these 
iV — 1 states do have dipole-allowed transitions to the 
second exciton band. They therefore have the poten¬ 
tial for preventing re-emission of light, while still allow¬ 
ing further absorption of photons, creating more excitons 


2 



FIG. 1. Left: N = 4 ring structure with attached trap. 
Right: Energy level diagram of this ring showing dipole- 
allowed optical transitions (red). Solid lines denote primary 
absorption transitions for the ratcheting cycle, others are in¬ 
dicated with thinner dashed lines. Only relevant phonon tran¬ 
sitions are shown (blue). Rarely populated levels are shown 
in grey. 

that in a photovoltaic circuit can be converted into useful 
photocurrent. Such states are examples of ratchet states 
\K r ), which have the simultaneous properties: 

t k r '■= K R ,K'5\K R \,\K'\+\ = °, (2) 

K' 

T K r : ='^2 T K r ,K'S\ Kr |,|a:'|-1 > °, (3) 

K' 

where is the sum of all transition rates connecting 

the ratchet eigenstate \Kr) to the adjacent band below 
(above) and Sij is the Kronecker symbol. 

For N = 4 and S > 0 the bright state in the single ex- 
citon band has the highest energy u + 25, and the three 
ratchet states have energies ui, uj and ui — 25, see Fig. 
and [IS] . In a molecular system, the ratchet states can be 
reached following excitation into the bright state when 
the exciton loses a small amount of energy to molecu¬ 
lar phonons. Importantly, such vibrational modes have 
much shorter wavelengths than optical modes and can of¬ 
ten be assumed to be local to each site. Unlike the photon 
field, the phonon interaction then breaks the symmetry 
and enables intraband transitions. 

Dynamics - In order to model the dynamics we use 
the QuTiP package [23] to construct full open system 
dynamics in Bloch-Redfield theory, for each kind of envi¬ 
ronmental interaction. We always retain the full form of 
the resulting dissipator tensor since the secular approx¬ 
imation cannot be applied when a system has bands of 
closely spaced levels which may be more closely spaced 
than the corresponding dissipation rates gUMj. 

First, consider the dynamics due to the light-matter 
interaction Hamiltonian 

N 

Hs-i = E E g q(^q+<E4) - ( 4 ) 

Q i =1 


where 6 q and Sq are creation and annihilation operators 
for photons of wavevector Q. The photon spectral den¬ 
sity is constant across all transition energies, and enters 
into the master equation dissipator only via the single 
emitter decay rate j a . Absorption and stimulated emis¬ 
sion terms are weighted by the appropriate Bose-Einstein 
factor N(un,T 0 ) = l/(exp[u>n/(k B T 0 )\ - 1), represent¬ 
ing the thermal occupancy of the optical mode with in¬ 
terband transition frequency wq and at temperature T 0 
(= 5800 K for solar radiation). 

Second, the exciton-phonon interaction is accounted 
for through [25] 

N 

H S -p = 'y ' y ' 9q,i&i (®q ,i + ®q,i) ) (5) 

2—1 q 

where for each site i a q and are the creation and an¬ 
nihilation operators for phonons of wavevector q with 
energy w q and exciton-phonon coupling g q . This inter¬ 
action commutes with of and so only mediates tran¬ 
sitions within each excited band. For simplicity, we again 
assume a spectral density J(w) that is approximately 
constant across all intraband transition energies wn, such 
that J(wn) = 7 p . We calculate appropriate phonon ma¬ 
trix elements for each transition, letting absorption and 
stimulated emission terms carry the appropriate Bose- 
Einstein factors N(u>u,T p ), with T p the ambient phonon 
temperature. Reflecting the fact that phonon relaxation 
typically proceeds orders of magnitude faster than optical 
transition rates [2D], we fix = 10007 o . This separation 
of timescales allows population to leave the bright states 
before reemission occurs. 

To complete the photovoltaic circuit, we include an 
additional ‘trap site’ f, with excited state |a) and ground 
state |/3), where excitons are irreversibly converted into 
work. The trap Hamiltonian is H t = w f af~of" , where 
uit represents its transition frequency, and of" = |/3) (a|. 
Henceforth, we shall consider the joint density matrix of 
ring system and trap, p = p a ® pt- Throughout the paper 
we use an incoherent hopping from a single site in the 
ring to the trap (see SI :I5J for an alternative model), 
and we fix the trap energy to be resonant with the state 
in which population is most likely to accumulate — unless 
otherwise stated, this is at the bottom of the first exciton 
band, Wt = w — 25. 

The theory of quantum heat engines [3 [23 [22] provides 
a way for us to assess the current and voltage output of 
the ring device: The action of a load across the device 
is mimicked by the trap decay rate 7 1 , which varies de¬ 
pending on the load resistance. The current is simply 
I = ey t(Pa)sa, with ( p a ) ss the steady state population of 
the excited trap state. The potential difference seen by 
the load is given by the deviation of the trap’s population 
from its thermal distribution [3121]: 

eU = ^ t + fc s T p ln^EE^ , 


( 6 ) 
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Parameter 

Symbol 

Default value 

Atomic transition frequency 

U) 

1.8 eV - 2S 

Hopping strength 

s 

0.02 eV 

Spontaneous emission rate 

7o 

1.52 x 10 9 Hz = 1 peV 

Phonon relaxation rate 

7p 

10 3 7o 

Extraction rate 

7* 

10 "So 

Photon bath temperature 

To 

5800 K 

Phonon bath temperature 

T P 

300 K 


TABLE I. Default model parameters; italic rows are plot pa¬ 
rameters in certain Figures, as stated in the relevant captions. 

a) b) 



FIG. 2. (a) Steady-state exciton population and (b) relative 

ratchet enhancement compared to the FD scenario (see text), 
both as functions of the optical and phononic bath temper¬ 
atures, and without any trap. The former is given in solar 
temperature units, K$ = 5800 K, other parameters are as in 
Table [T] except that 7 ^ = 0. 

where ks is Boltzmann’s constant and T p the (ambient) 
phonon bath temperature. To study the relationship be¬ 
tween current, voltage, and therefore power, we alter the 
trapping rate 7 *. Letting 7 t —> 0 or 00 leads to the well- 
known open and short circuit limit, respectively, however, 
our main interest in the following is in the region near 
the optimal power output point of the device. 

Results - We proceed by finding the steady state of the 
following master equation which incorporates the optical 
and phononic dissipator as well as the trap decay process: 

p = —i[Hs + Ht, p\ + D 0 [p] + D p [p] + D t [p] + D x [p \, 

(7) 

where D„[p ], D p [p\ are the Bloch Redfield dissipators 
for the photon and phonon fields discussed above, and 
D t [p\ = 7 t(<r t _ paf — ^{cr^eqr, p}) is a standard Lindblad 
dissipator representing trap decay. D x [p\ represents ex¬ 
citon extraction via incoherent hopping at rate "f x from 
a single site on the ring to the trap (see Fig. [I] a). 

To expose the effect of the optical ratchets we now 
contrast the full model of Eq. 0 with two other artifi¬ 
cial cases that exclude the possibility of ratchets: In the 


first, we remove all phonon-assisted relaxation (y p = 0 ), 
meaning the system only undergoes optical transitions 
along the Dicke ladder of bright states [4] - we term this 
‘no phonons’ (NP). In the second construct, the ratchet 
states are rendered fully dark by setting all upwards tran¬ 
sition matrix elements from them to zero, a scenario we 
refer to as ‘forced dark’ (FD). All results in this section 
are based on the steady state of the system for N = 4, ob¬ 
tained by an iterative numerical method performed with 
QuTiP |23j . and using the default parameters from Ta¬ 
ble |T] unless explicitly stated otherwise. Our choice of 
uj = 1.8 eV—2 S differs from that in previous works 
our subtraction of 2 S ensures that the highest (bright 
state) level in the first band always remains at 1.8 eV 
and absorbs at a fixed frequency. 

In Fig.[2]we display the steady-state exciton population 
of the system without trapping as a function of optical 
and phonon bath temperatures. We compare the ratchet 
model to FD, finding that the ratchet enhancement is 
greatest in the hot photons, cold phonons regime. This is 
to be expected: hot photons quickly promote the system 
up the excited bands via ratchets, with cold phonons al¬ 
lowing a one-way protection of gained excitation energy. 
For FD states at low phonon temperature the exciton 
population plateaus near one, as most population ends 
up trapped in the state at the bottom of the first band 
(see HU for the FD data and an extended discussion). 
By contrast, the ratchet states keep absorbing, allowing 
the steady state population to rise toward the infinite 
temperature limit of 2 (for N = 4). 

Using the heat engine model, we can now explore the 
performance of the four site ring as an energy harvester. 
We will focus on the dependence of the power output on 
two parameters: the coupling S, and extraction rate 7 X . 
For each parameter set, the trap decay rate 7 1 is varied 
for maximum power output; example power traces as a 
function of 7 1 can be found the SI CU- 

In Fig. [3] we display the absolute value of the power 
generated for the ratchet and NP scenarios, whereas in 
Fig. a we display the relative power output of ratchet 
states over each of the two artificial scenarios. To en¬ 
sure fair comparison, for NP we set the trap energy to 
be resonant with the bright state at w + 2 S, since the 
ratchet states at lower energy are inaccessible here. All 
three scenarios have lower output power when the trap¬ 
ping (charge separation) rates are low, since this creates 
a bottleneck in the cycle and limits the size of the pho- 
tocurrent. 

However, the ratchet states hold a great advantage over 
the other scenarios in this bottlenecked region, since this 
is precisely the situation in which excitations need to 
be held for some time by the ring before extraction is 
possible. The NP scenario only performs poorly in this 
case, since excitations in all likelihood decay before being 
extracted. Indeed, the relative power output of ratchet¬ 
ing over NP, shown in the left panel of Fig. [4j therefore 
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FIG. 3. Absolute power output in the ratcheting (blue sur¬ 
face) and NP (red surface) cases as a function of (single site) 
extraction rate 7 ^ and hopping strength S. Other parameters 
are as in Table I. For each point, the optimal trap decay rate 
is found by numerical search. The sets of curves projected 
onto the side walls of the 3D plot are cuts through the data 
for fixed values of the appropriate parameter. The two sur¬ 
faces cross in the region where the bottleneck is lifted, with 
ratcheting generating more power for a severe and intermedi¬ 
ate bottleneck. The line at which the curves cross is projected 
onto the bottom (xy) plane of the figure. 

rises to as high as a factor of 35. The right panel of 
f^. a demonstrates the importance of ratcheting over 
dark state protection alone, with the ratchet model de¬ 
livering up to 20 % better performance than FD. 

In Fig. a w e can see that advantage of ratcheting per¬ 
sists into a region where there is only a moderate bot¬ 
tleneck, i.e. up to ~ 10 -5 eV = IO 70 . Referring to 
Fig. awe see that the ratchet power output in this region 
is already close to its maximum. Even so, if the extrac¬ 
tion rate was arbitrarily tunable, then even a moderate 
bottleneck could be avoided completely and the main ad¬ 
vantage of ratcheting removed. However, in photosynthe¬ 
sis, creating a fast extraction rate carries with it a clear 
resource cost: whereas antenna systems can be compar¬ 
atively ‘cheap’, reaction centres carry a much larger spa¬ 
tial footprint, being typically embedded in membranes 
and requiring significant surrounding infrastructure (such 
as concentration gradients produced by proton pumps). 
The severity of the bottleneck is likely to be inversely pro¬ 
portional to the number of reaction centres. A photocell 
design exploiting ratcheting in the moderate bottleneck 
regime would be likely then to generate optimal power 
per unit volume of material — and similar design prin¬ 
ciples will apply to artificially designed molecular light 
harvesting systems. 

The choice of interaction strength S is also important. 
Ratcheting achieves optimal results for S ~ 0.05 eV in 
the bottlenecked region. This dependence arises because 
the size of S determines the gap between bright and dark 
or ratchet states. In turn, this controls the effective rate 



FIG. 4. Enhancement of power output for ratchet states over 
NP (left, displayed as a multiplicative factor) and FD (right, 
displayed as a percentage), as a function of both hopping 
strength S and extraction rate -y x . Parameters are as in Fig. [3] 

for ‘upwards’ phonon-assisted transitions within excita¬ 
tion bands, as those rely on the absorption of a phonon 
and are proportional to N(lou,T p ). Consequently, larger 
values of S entail more directed dissipation into the lower 
states of each band, boosting the occupation of ratchet 
states. However at the same time, increasing S leads to 
a lower ratchet state energy and so a lower trap energy - 
and hence to a voltage drop. The trade off between these 
two competing influences leads to a maximum in ratchet 
performance. 

As we discuss in detail in the SI, ratcheting continues to 
convey an advantage in the presence of moderate levels of 
various real-world imperfections, such as site energy dis¬ 
order, non-radiative recombination, and exciton-exciton 
annihilation m- 

Conclusion - We have investigated the light harvest¬ 
ing properties of coupled ring structures, inspired by the 
molecular rings that serve as antennae in photosynthesis. 
Considering a vibrational as well as an electromagnetic 
environment allows the system to explore the full Hilbert 
space rather than just the restricted subset of Dicke lad¬ 
der states. We have shown that the off-ladder states pos¬ 
sess interesting and desirable properties, which can be 
harnessed for enhancing both the current and power of a 
ring-based photocell device. 

Several possible systems m could be used to observe 
the effect, from superconducting qubits m to macro- 
cyclic molecules [23]. Our approach generalises exist¬ 
ing concepts for dark state protection to arbitrary 
numbers of sites and importantly includes multi-exciton 
states, which introduce the ratcheting effect as a distinct 
additional mechanism for enhancing the overall light¬ 
harvesting performance. The optical ratchet enhance¬ 
ment is particularly well suited to situations where ex- 
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citon extraction and conversion represent the bottleneck 
of a photocell cycle. 

In future work it would be interesting to explore com¬ 
bining optical ratcheting with enhancements of the pri¬ 
mary absorption process, for example by exploiting the 
phenomena of stimulated absorption [28j or superabsorp¬ 
tion g]. 
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Supplemental Information: 
Quantum-enhanced capture of photons 
using optical ratchet states 


RING STRUCTURE DIAGONALISATION 


The Jordan-Wigner transformation maps a Pauli spin 
1/2 system onto a ‘hard-core’ boson model [ST]. This 
leads to a bosonic description of the collective excitons 
whilst maintaining Pauli’s exclusion principle, which for¬ 
bids double excitation of a single site. In our case each 
‘spin’ represents one of N identical optical emitters / 
absorbers and the spin’s ‘up’ / ‘down’ projection along 
the 2 -axis denotes the presence / absence of an exciton 
on the respective site. We assume a ring-like geometri¬ 
cal arrangement as shown in Fig. la of the main paper. 
Considering a ring rather than a spin chain with free ends 
introduces the additional complexity of an alternating pe¬ 
riodic or anti-periodic boundary condition depending on 
the number of excitations jS2] , 

The Jordan-Wigner transformation is achieved by in¬ 
troducing bosonic creation and annihilation operators, 
which are defined in terms of the Pauli spin raising and 
lowering operators as follows: 


c[ =a + e ^r =1 a. 


Cl 


E l — 1 

_ 3=1 


J t a i <j7 


(51) 

(52) 


Here, the e I7r2jl factors are traditionally referred 

to as ‘strings’, and their inclusion is necessary for produc¬ 
ing the correct anti-commutation relation for fermions, 
which also apply to hard-core bosons: 


{c a , cj)}— ; (S3) 

{ 4 : 4 }= {c«, C^} = 0. (S4) 


Performing the transform (Si] S2) to the Hamiltonian of 
the main paper given by Eq. (1) yields: 


N 


H s =wh- S J2( c ]cj+i + c] +1 Cj) 


i= i 


+ S(e l7Tn + l)(cjyCi + c-^cn) . 


(S5) 


where h = J2f= l c ) c i ~ 


_ N 

- 2^7 = 1 


u1 (jj counts the num¬ 


ber of excitons on the ring. The final term in Eq. (S5) 


has been introduced to implement the required alternat¬ 
ing boundary condition [32]. Note that this elegant and 
compact form of the Hamiltonian double-counts the in¬ 
teraction term for N = 2, but this is of no concern here 
since we are only interested in N > 4. 

Cle arly , the operator e I7m commutes with the Hamilto¬ 
nian (S5|, so that the eigenvalues of e I7m are good quan¬ 


tum numbers. The system can thus be partitioned into 
two parity subspaces of even and odd exciton number n , 


each of which can be solved separately. We complete the 
diagonalisation with the help of the Fourier transform 


° k — Vn e C U 

(S6) 

r , — 1 r 1 p-ikjr 

Ck - vW 2-b=i e ’ 

(S7) 

where the kj are phase factors defined by 


k j = W3 ( n odd ) > 

(S8) 

kj = ( n even ) ■ 

(S9) 


Each of the above N distinct fc-values corresponds to 
a single exciton state of the system. The eigenstates of 
the full multi-exciton system are given by combinations 
of these states. The eigenstates are thus defined by sets 
K whose elements are the fc’s, and where the number of 
elements in K represents the number of excitons in the 
state ( \K\ = n). The Pauli exclusion principle manifests 
itself by allowing only different single exciton states to 
form a multiple exciton state. The eigenvalues are then 
given by [ST 


A k = X (w — 2 S cos fc) 

k£K 

with corresponding eigenstates 
- J 


(S10) 


1 |0>, (S11) 


1*0 = c * !1 - c fc. 

where |0) is the ground state. Obviously, Eq. ( Sll| ) yields 
an eigenstate with n excitations and each band comprises 
a total of n! (4-rt)! P oss ible states. In its diagonalised 
form the Hamiltonian Eq. ( |S5| ) is then simply given by 

H s = J2^k \K)( K \. (S12) 

K 


OPTICAL TRANSITON MATRIX ELEMENTS 

The optical transition operator ,4 = X)4i <4 con¬ 
nects eigenstates differing by one exciton. The transition 
rates are proportional to T KK / = \{K'\J + \K)\ 2 and can 
be explicitly evaluated with the help of Slater determi¬ 
nants, albeit involving some tedious calculations fS2[ [S3], 
yielding: 


k,k' — 


2«aK-"+5)j ( 


(S13) 


IL i >i'(e ik *-e* k *')n j>j ,(e- ik i-e~ ik i') 






) 


where S(x,y) is 1 if x = y + 27r to for integer m and 
zero otherwise. The indices inside the sums and 

products of the above expression refer to the elements of 
unique pairs of fc-values composing I\ or K' (see |S2| for 
full details). 
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Ratchets steady-state number of excitons 


FD steady-state number of excitons 


NP steady-state number of excitons 



FIG. SI. Number of excitons in the ring molecule in the steady state, with no attached trap, (a) Ratchets model; (b) forced 
dark (FD) model; (c) no phonon (NP) model. All three panels use the same colour mapping for ease of direct visual comparison. 
Panel (a) plots the same data also displayed in Fig. 2a in the main text. Parameters are from Table I in the main text. 


STEADY STATE EXCITON POPULATIONS 


In Fig. 2 of the main text we display the steady-state 
exciton population in the ratchets model as a function 
of phonon and photon temperatures, and also the rel¬ 
ative population of the ratchets model over the forced 


dark (FD) scenario. To complete the picture, in Fig. SI 


we show the absolute values of the steady-state exciton 
count for all three scenarios. In the no phonons (NP) 
case, as expected there is no dependence on the phonon 
temperature. Indeed, since there is no dark-state protec¬ 
tion in the NP case, we find a lower exciton population 
than in the other two cases; the difference is more pro¬ 
nounced at lower phonon bath temperature, when dark 
states are most likely to be populated. Naively, one might 
expect the surface for the FD scenario to plateau at one, 
however, the ladder of bright states remains available 
in this case and via it population can access and get 
trapped in higher lying dark states above the single ex¬ 
citon band. We expect this to happen for high optical 
and low but finite phonon temperature, and indeed in 
this region the surface rises above one. Note that unlike 
panels (a) and (c) in Fig. [STJ (b) starts at a phonon bath 
temperature of T = 75 K due to numerical instabilities in 
our steady-state solver. At much lower phonon temper¬ 
ature, where phonon-assisted relaxation becomes effec¬ 
tively uni-directional, we expect the FD surface to drop 
towards one, as it gets increasingly unlikely for popula¬ 
tion to bypass or escape the single excitation subspace. 
Using alternative methods, we have verified this is indeed 
the case (not shown). 


0.15 



ng Sl~r~ 0.03 


FIG. S2. Absolute power output in the ratcheting (blue sur¬ 
face) and FD (red surface) cases as a function of (single site) 
extraction rate and Forster coupling. Here the collective ex¬ 
traction model is used. Other parameters are from Table I 
and thus identical to those of Fig. 4 of the main text. The 
sets of curves projected onto the side walls of the 3D plot 
are cuts through the data for fixed values of the appropri¬ 
ate parameter. For each point, the optimal trap decay rate 
is found by numerical search. The two surfaces cross in the 
region where the bottleneck is lifted, with ratcheting generat¬ 
ing more power for a severe and intermediate bottleneck. The 
line at which the curves cross is projected onto the bottom 
(xy) plane of the figure. 
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FIG. S3. For the collective extraction model, enhancement 
of power output for ratchet states over NP (left, displayed as 
a multiplicative factor) and FD (right, displayed as a percent¬ 
age), as a function of both hopping strength 5 and extraction 
rate The simulations use the same model and parameters 
as Fig. [S2] 

EXTRACTION MECHANISMS 

Our photocell cycle relies on the extraction of the ex- 
citons from the ring to the trap site. Below we give an 
explicit description of two different mechanisms and mod¬ 
els, one of which features in the main text, the other 
whose results are displayed for comparison in this sup¬ 
plement. 

For the incoherent hopping process that we use in the 
main text we assume that excitons are extracted incoher¬ 
ently from a single site of the ring to the trap. We then 
simply include a dissipator: 

D x [p] = cr+, p\ . (S14) 

For optimal extraction the trap’s energy is resonant with 
the desired extraction level, i.e. for the ratcheting and 
FD cases, the lowest state in the first exciton band, w t = 
uj—2S] for the NP case, the trap is instead made resonant 
with the bright state, u>t = ui + 2 S. For consistency and 
a fair comparison across the different mechanisms, the 
energy of the trap site is also adjusted in the same way 
for the extraction mechanism described next. 

We use the single site extraction process in the pa¬ 
per since it is relatively straightforward to implement. 
However, it is not necessarily the most efficient — so 
in this supplement we will also consider the possibility 
of extraction directly from particular eigenstates of the 
ring. This can sometimes be achieved by suitable posi¬ 
tioning of the trap: e.g. placing the trap at the ring’s 
centre with equal coupling to all other sites allows ex¬ 
traction from symmetric states |S4| . Unfortunately, this 
does not straightforwardly generalise to the more com¬ 


7r ( Hz ) 

10 16 10 14 10 12 10 10 10 s 10 6 10 4 10 2 



FIG. S4. Power (solid) and its limit (dashed) as a function 
of voltage based on single site incoherent Lindblad extrac¬ 
tion. Three scenarios are shown, ratcheting (blue), forced- 
dark (FD, green) and NP (red). The power limit is given 
by the maximum possible voltage multiplied by the current 
P max = I(u> + 25) IS5I . Model parameters are from Table I 
of the main text. Each plot is found by varying yt indepen¬ 
dently for the different models. The upper *-axis shows how 
7 1 varies in the ratchets model over many orders of magnitude 
for relatively small changes on the voltage scale. 

plex anti-symmetric dark states. For the special cases 
of N = 2, 3 alternating phases on the interaction terms 
achieves the desired coupling [S5; jS 6 j; however, scaling 
up to more sites would require increasingly exotic molec¬ 
ular orbitals. We therefore implement a collective state 
extraction method using a set of operationally motivated 
Lindblad operators C). In the case of the ratchets and 
the FD model, the C\ move population from the lowest 
state of a particular band of H s into the lowest state in 
the exciton band below, exciting the trap in the process. 
In the NP case, the Ci go from the highest state of one 
band to the highest of the band below (since these are 
the only states accessed in this model). We incorporate 
this process into our model through the dissipator 

D x , a d\p]=’yx'52v[C i ,p], (S15) 

i 

where, once again, 7 X is the effective extraction rate. 
Once on the trap, excitons are assumed to undergo charge 
separation at rate 7 1 , which we model as an irreversible 
decay process. 

Fig. [S 2 | shows the absolute power for the ratchets and 
NP model for the collective trapping model. This plot is 
equivalent to Fig. 3 of the main paper which displayed 
the same quantities for the single site incoherent trap¬ 
ping model. Similarly, Fig. [S3] shows the relative power 
enhancement of ratcheting over both artificial models for 
the collective trapping model — again this is the equiv¬ 
alent plot to Fig. 4 in the main paper but here for the 
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FIG. S5. Power-Voltage characteristics for the single site incoherent extraction model. Parameters are as in Table I of the 
main paper except for 7 ^. Each panel compares the three models (ratchets, NP, FD) for a different extraction rate 7 ^. The 
advantage of ratcheting is most apparent for slower extraction rates, i.e. in the bottleneck region. 


collective extraction model. Unsurprisingly, the absolute 
power extracted in the collective model is larger than the 
incoherent model. This is because we are here extracting 
directly from the ratchet or dark states and thus have ac¬ 
cess to the ‘entire’ exciton rather than just its projection 
onto a single site of the ring. The relative advantage of 
ratcheting over FD has a similar form for both extraction 
mechanisms, though it is smaller for the collective model. 

POWER-VOLTAGE CURVES 

In Figs. 3 and 4 of the main paper we are displaying 
maximum output power as a function of S and "f x . Each 
point in those plots is a result of generating a power- 
voltage curve for each choice of parameter, and finding 
the maximum power. Each such curve is generated by 
varying the trap decay rate 7 1 - We show examples of 
these power-voltage curves in Fig. [S4j considering the 
usual three scenarios: ratchets, FD and NP. On the upper 
x-axis we also show the values of 7 * corresponding to each 
voltage taking the example of the ratchets case; as shown 


7 t varies over orders of magnitude for only modest voltage 
shifts. These curves make use of the parameters specified 
in Table I of the main text, and a comparison with Figs. 3 
and 4 of the main paper shows that S = 0.02 eV and 
7 a, = 10~' eV are in fact not optimal for showing the full 
ratcheting advantage. 


Further examples of P V curves are shown in Fig. |S5[ 
each for different single site incoherent extraction rates, 
fixing S = 0.02 eV. We can see here the advantage of 
ratcheting in the intermediate and severe bottleneck re¬ 
gions. For very fast extraction, any absorbed photon gets 
immediately transferred to the trap, so that neither dark- 
state protection nor ratcheting offer any advantage (the 
system is already operating as efficiently as possible). In 
this case NP achieves the same current as ratchets and 
FD, but ends up producing a higher resultant power as 
its trap transition frequency - and thus effective voltage 
- is higher, as explained in Sec. [J 
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FIG. S 6 . Left: absolute output power as a function of S 
(for fixed 7 ^ = 10 -7 eV). Right: absolute output power as a 
function of 7 X (for fixed S = 0.02 eV). Curves for all three 
scenarios - ratchets (R), NP, FD - are displayed. The two 
extraction mechanisms of incoherent single site (inc) and in¬ 
coherent collective (col) are also shown. All other parameters 
are as specified in Table I of the main text. 


EXTRACTION BOTTLENECK 2D PLOTS 

In Fig.[S 6 ]we show several 2D sections of the surfaces 
in Fig. 3 of the main paper and Fig. [S2] i.e. we show 
absolute output power as a function of S (for fixed 7 ^ = 
10 - 7 eV) in the left panel and as a function of 7 ^ (for 
fixed S = 0.02 eV) in the right panel. These are displayed 
for all three scenarios (ratchets, FD, NP) and for both 
collective and incoherent extraction mechanisms. 


IMPERFECTIONS 


Disorder 


To model disorder, we generate random site energies 
that are normally distributed around the transition fre¬ 
quency oj + 2S = 1.8 eV with standard deviation < 7 . The 
/, V value pairs for each trap decay frequency are av¬ 
eraged over 100 different implementations, and are based 
on the Bloch-Redfield dissipator [S7] for the optical and 
phonon interactions. 

Results for different degrees of disorder are shown in 


the left of Fig. [S7j For small levels of disorder er < S, 
there is little deviation from the ideal case, and 100 im¬ 
plementations are indeed enough to give fairly smooth 
curves. Once a approaches S (0.02 eV corresponds to a 
a of around 1.14%), however, not only is the achievable 
power with ratchets visibly reduced but our statistical 
ensemble is also not large enough to smooth out those 
larger performance variations, resulting in jagged-looking 
curves. 

Interestingly, the performance of the NP case improves 
with increasing disorder. This may be due to the fact that 
even in asymmetric circumstances eigenstates with a var¬ 


ied range of interesting optical properties may arise jS8] 
in the fully accessible Hilbert space, and the absorption 
dynamics is no longer restricted to only the single lowest 
ladder transition. 

We do not show the FD scenario here, as the pres¬ 
ence of disorder makes it impossible to properly distin¬ 
guish between allowed ‘ladder’ and FD ‘forbidden’ opti¬ 
cal transitions. However, we note that random sampling 
for small levels of disorder (well below relative site en¬ 
ergy fluctuations on the order of a percent) - where the 
symmetries and properties of system eigenstates remain 
largely unperturbed - indicates that FD always performs 
worse than the full ratchets model, as expected. 

The right panel of Fig. [S7] shows the averaging from 
a single up to 1000 implementations of randomness for 
a = S/2. Interestingly, the single run has some data- 
points which can even exceed the ideal, fully ordered, 
ratcheting performance — though the ensemble average 
power is clearly below that in the fully ordered case. We 
suspect these outliers correspond to highly efficient asym¬ 
metric configurations: For the case of a dimer and dark- 
state protection it has recently been established that cer¬ 
tain asymmetric arrangements can potentially outper¬ 
form symmetric ones [S8l . However, the unravelling of 
the precise interplay of mechanisms at play will be more 
complicated in this larger system, but would be an inter¬ 
esting line of investigation for a future study. 

Note that we keep our extraction and trap model unal¬ 
tered for our disorder calculations, i.e. we extract with a 
trap energy of w t = 1.8 eV ± 2S for the NP and ratchets 
case, respectively. This simplifying choice may in general 
entail slightly over- or underestimating the effective volt¬ 
age of our trap. However, the relative voltage discrepancy 
arising from this assumption is expected to be limited by 
energy shifts due to disorder given by a. As seen in both 
panels of Fig. [S7j disorder in the model affects the power 
output for both ratchets and NP much more significantly 
than that, meaning our results are clearly dominated by 
other processes as opposed to just being an artefact of 
our simplified trapping model. 


Non-radiative decay 

We model non-radiative decay by including a full set 
of uni-directional Lindblad lowering processes acting in 
the site basis (as opposed to the collective optical decay 
events which are present due to the interaction with the 
shared electromagnetic environment). 

Increasing non-radiative decay rates significantly de¬ 
creases the achievable current and power output of 
the photocell as shown in Fig. |S8| Interestingly, non- 
radiative decay does not affect all three scenarios equally, 
and the ratcheting advantage diminishes with larger non- 
radiative recombination rate. The ratchets states enjoy 
no special protection against non-radiative decay act- 
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FIG. S7. Effect of site energy disorder on the achievable power. Left: The darkest blue and red curves represent the ideal case 
of no disorder, the lighter curves with dots correspond to random site energy fluctuations with a standard deviation from the 
idealised case of 0.1%, 0.33%, 1%, and 2%. All these curves are averaged over 100 different randomly selected configurations. 
Right: Convergence of disorder curves for a = S/2 (corresponding to about 0.57% disorder of site energies). The solid lines 
represent the average over 1000 runs, whereas the joined circles from faint to more opaque correspond to the average of 1, 10, 
and 100 runs. Note that certain random configurations prove to be highly efficient - see text for a discussion. Model parameters 
are from Table I of the main text. 
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FIG. S8. Non-radiative decay processes acting independently on each site, with effective rates as shown above each panel. All 
other parameters are from Table I of the main text. Aggressive non-radiative recombination severely reduces the achievable 
power output of all three models, ratchets, NP, and FD, lending an edge to NP where excitations are never stored for long in 
excited levels. 


ing on individual sites, and when they carry more than 
one excitation the impact of non-radiative decay is felt 
more acutely. However, the combined ratcheting and 
dark state protection advantage persists as long as non- 
radiative decay does not become the dominant loss chan¬ 
nel. 


Exciton-exciton annihilation 

We now discuss the impact of exciton-exciton annihi¬ 
lation (EEA) on the ratcheting efficiency enhancement. 
EEA is the result of a two step dynamical process [55] , 
First, excitons in adjacent sites on the ring can fuse to¬ 
gether, creating an excitation on a single chromophore 
that has twice the energy of the first excited exciton state 
so far considered. The mechanism for this is, like exciton 
transfer, Forster coupling - i.e. a dynamic dipole-dipole 
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FIG. S9. The effect of exciton-exciton annihilation at an effective rate as shown above each panel. All other parameters are 
from Table I of the main text. As expected neither NP nor FD are much affected, and notably the ratcheting advantage only 
fully vanishes once the exciton-exciton annihilation rate exceeds the spontaneous emission rate by an order of magnitude or 
more. 


interaction between the adjacent excitations. The sec¬ 
ond step is fast non-radiative decay of the double energy 
state on the excited chromophore back down to the reg¬ 
ular single exciton state. 

The Hamiltonian underlying the first step might be 
expressed by using an extension of the notation we have 
employed in the main paper. We now allow three levels 
per chromophore: the already-considered ground state 
|G) and low energy exciton state | E), as well as the 
newly introduced double energy excitation | D). With 
these state labels the Pauli operator we use in the main 
text is cr + = \E) (G|, and now we define a new spin raising 
operator rj + = \D) (E\. The Hamiltonian then becomes: 

N 

H s' =u °i + 2 ViVi )+ (S16) 

i—1 

N 

s H( CT * +cr i+i + a tvT+1 + H c -) > ( S17 ) 

i—l 


where we have made the simplifying assumption that the 
two excitation energy transfer has the same strength (S) 
as the single excitation one. This is reasonable, implying 
a similar transition dipole from |G) to | E) as from | E) 
to |G). By diagonalizing Eq. (S16) for an iV = 4 system, 
in a basis of states {|G), | E ), | D)} on each site, we are 
able to straightforwardly obtain the new eigenstates of 
the system. Again, since the Hamiltonian preserves to¬ 
tal excitation number (assuming that | D) states count as 
two excitations), then the eigenstates fall into a series of 
bands, each corresponding to a different excitation num¬ 
ber. The zero excitation ground state and the four single 
excitation eigenstates are obviously unaffected by the | D) 
states, but the second band now consists of ten, not six, 
possible states, which are now combinations of both | D) 
states on single sites and | E) states on two sites. The 


lowest energy state in this second band (let us label it 
| Q) 2 ) would be reached following rapid relaxation follow¬ 
ing ratchet action from the first band. We have calcu¬ 
lated that | Q) 2 contains only a 16% |D) state character 
(by probability). Considering that population of the | D) 
state is crucial for EEA to occur, we can conclude that 
EEA will be significantly suppressed from | Q) 0 . 

Moreover, we have calculated the electric transition 
dipole moment between | Q) 2 and all states in the first 
band. The dipole is zero for transitions into all but one 
of the first band states, and only 18% of the dipole of 
a single uncoupled chromophore for the remaining state. 
By contrast, if the | D) states are not included, the equiv¬ 
alent transition dipole is 41% of the single chromophore 
value. Thus, we can conclude that the lowest state of the 
second band is more protected from optical re-emission 
if the | D) states are included in the model. This means 
that including the possibility of doubly-excited sites sur¬ 
prisingly lends a certain degree of extra ‘dark-state pro¬ 
tection’ character to the second excitation band, allowing 
more time for exciton harvesting, and thus improving the 
efficiency of the device. 

There are materials that do not quickly decay through 
non-radiative transition to the lowest energy state of a 
given band - these materials then do not obey Kasha’s 
rule jSlOj . which states that optical emission is always 
from the lowest energy state of a band. For these ma¬ 
terials, the reasoning above regarding the composition 
of the lowest second band state is less relevant. How¬ 
ever, in that case, the second step involved in EEA - the 
non-radiative decay of the | D) state - is also inhibited, 
and this too leads to an inhibition of EEA and so to 
more efficient energy harvesting devices. Indeed, finding 
materials that violate Kasha’s rule for novel optoelec¬ 
tronic applications such as this is a very active area of 





























research mm- 

As discussed above, we expect the | Q) 2 to have sup¬ 
pressed EEA, but we do not expect it to disappear com¬ 
pletely. We therefore look at its impact by introducing 
EEA it into our model by means of one-way Lindblad 
operators taking population from all doubly or higher 
excited states to the highest energy (ladder) state in the 
single excitation subspace. In Fig. |S9| we show power- 
voltage characteristics with such EEA taken into account. 
Each plot compares the ratchets, FD and NP models, for 
the value of the annihilation rate indicated above it. 

As expected, neither the NP case nor the FD scenario 
are visibly affected. As annihilation increases, the maxi¬ 
mal current achieved by the ratchets case eventually re¬ 
duces to that of FD but, importantly, the ratcheting 
advantage only fully vanishes once the exciton-exciton 
annihilation rate exceeds the spontaneous emission rate 
by an order of magnitude or more. On close numerical 
inspection, FD and NP also suffer very slightly under 
EEA, because in both those cases the ladder of bright 
states remains available and there is thus slight chance 
of multiple excitation. However, unlike in the case of 
non-radiative recombination, annihilation never kills the 
dark-state protection advantage. 

DEMONSTRATION SYSTEMS 

In this section we mention some possible systems in 
which the action of ratchet states might be initially de¬ 
tected. 

Ideally, we require a system composed of at least four 
coupled, radiation-absorbing elements, and where the 
spacing between each element is either much less than 
the wavelength of light, or else where the elements can 
be placed such that the radiation with which they inter¬ 
act is in phase with all of them. This latter condition 
allows for the possibility of using systems placed an in¬ 
teger number of wavelengths apart, and would include, 
for example, superconducting qubits interacting with mi¬ 
crowave cavity fields. For example, two such qubits have 
been used to demonstrate super-radiance |S13j . Includ¬ 
ing two more would enable ratchets states to be reached, 
initially by driving the bright excitation and then by al¬ 
tering single qubit phases through local manipulation. 

A system closer in spirit to what we have described 
here would involve coupled optical absorbers. For exam¬ 
ple, NV _ centers in diamond are established as systems 
on which to test quantum information protocols |S14j 
and fundamental features of quantum theory [S15j . NVs 
can be individually addressed optically, and have spin se¬ 
lective coupling to light that allows optical measurement 
of their spin state. Moreover, they have been produced 
in patterns with a precision of around A/10, where A is 
the wavelength at which they emit Sl(i . The prospect 
of collective photon emission (superrradiance) with mul¬ 


tiple NV centres has recently been studied m- 

NVs have strain dependent splitting so would have 
some disorder, though as we have seen in Section | above, 
this does not mean that ratchets states lose all function. 
However, it may be that silicon vacancies (SiVs) in dia¬ 
mond ism. or alternatively (di)vacancies in SiC, make 
better test systems. For example, SiVs in diamond can 
be placed with A/20 precision and are very stable, with 
inhomogeneous broadening heavily suppressed by their 
inversion symmetry. For all these kinds of emitter, the 
Forster coupling could be enhanced using a waveguide or 
cavity to promote the interaction of several centers with a 
common optical mode |S19j . Transient luminescence and 
absorption techniques could be used to track the light 
excitation in time as it moves through the system, and 
demonstrate behavior consistent with ratcheting. 

Quantum dots are another system that can ex¬ 
hibit strong optical coupling, and could for example 
be arranged in sub-wavelength patterns using DNA 
origami jS20l . Evidence of collective photon emission of 
lateral quantum dot ensembles has also already been re¬ 
ported [S2l| . 

Of course, the final device likely to be used in a real 
photocell would be a molecular system. Such systems 
can be made with high reproducibility, and such that 
they are very closely spaced with respect to their emis¬ 
sion wavelength. H-aggregates of chromophores have a 
positive Forster coupling that shifts the bright (absorb¬ 
ing) state up in energy, as is the case in our model [S22] . 
Alternatively, light-harvesting 7r-conjugated supramolec- 
ular rings can be made with twelve or more repeating 
units |S23| , specifically to mimic the ring structures found 
in biology. Indeed, the symmetry of this structure gener¬ 
ates a dark Sq—Si transition in the absence of vibrational 
distortion [S24j . Further design optimization of both the 
individual systems and rigidity of the ring could result in 
a system ready to meet the requirements for ratcheting. 
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